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THE  STRUCTURE  OF  SHOCK  FRONTS  IN  ARGON  AND  NITROGEN* 


by 

Donald  F,  Homig  and  Melvin  Linzer'*’ 

Frick  Chemical  Laboratory 
Princeton  University,  Princeton,  N,  J, 

Abstract 

The  optical  reflectivity  method  was  used  to  investigate  the  structure 
of  shock  fronts  in  argon  from  Mach  1.70  to  Mach  4.85  and  in  nitrogen  from 
Mach  2.01  to  Mach  3«72.  Experimental  data  were  obtained  at  two  wavelengths 
and  over  a  wide  range  of  initial  pressures.  The  reflectivities,  corrected 
empirically  for  shock  curvature,  were  fitted  to  a  bimodal  profile  to  yield 
a  maximum«slope  density  thickness.  The  reciprocal  of  the  thickness  in  argon 
(expressed  in  terms  of  the  Maxwellian  mean  free  path  in  the  undisturbed  gas) 
rises  rapidly  to  a  maximum  of  approximately  0.31  at  about  Mach  3*5  and 
decreases  gradually  thereafter.  Above  Mach  3  the  thickness  is  about  5^^ 
greater  than  calculated  from  the  Navier-Stokes  equations,  using  a  realistic 
viscosity-temperature  relationship.  There  is  excellent  agreement,  especially 
at  the  higher  shock  strengths,  with  recent  bimodal  calculations  carried  out 
by  Muckenfuss,  using  realistic  intermolecular  potentials.  In  nitrogen,  the 
shocks  are  thinner  than  in  argon  and  appear  to  attain  a  minimum  value  of 
2.3  initial  mean  free  paths  at  about  Mach  3*7*  Rotational  relaxation  appears 
to  be  as  rapid  in  the  strong  shocks  as  previously  observed  in  weak  shocks; 
it  appears  to  be  completed  within  the  shock  front.  The  experimental  density 
thicknesses  are  approximately  5^^  greater  than  those  calculated  from  the 
Navier-Stokes  equations,  using  experimental  shear  viscosity  and  a  bulk 
viscosity  of  2|t/3.  The  agreement  with  these  Navier-Stokes  solutions  is 
about  as  good  as  those  in  argon. 

*The  work  reported  here  is  based  on  a  thesis  submitted  by  Melvin  Linscr  in 
partial  fulfillment  of  the  requirements  for  the  degree  of  Doctor  of  Philosophy 
at  Princeton  University,  1962. 

^Present  addresst  National  Bureau  of  Standards,  Washington,  D.  C. 


Introdciotiofi 


The  straeture  of  shock  fronts  in  various  gases  has  been  investigated 
by  a  number  of  experimental  techniques  in  the  last  few  years^*^*.  Up  to 
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about  Mach  2  no  clear  cut  differences  have  appeared  between  the  shock 
thicknesses  determined  experimentally  and  those  predicted  by  the  Navier. 
Stokes  equations,  a  surprising  result  in  view  of  the  large  deviations  from 
thermal  equilibrium  which  exist  within  the  shock.  Indeed,  it  has  been 
pointed  out  by  Sherman  and  Talbot^^  that  the  kinetic  theory  which  leads 
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to  the  Navier.Stokes  equations  breaks  down  completely  above  Mach  5*  giving 
a  physically  impossible  molecular  velocity  distribution. 

The  effect  of  higher  order  terms  in  the  kinetic  theory  has  been  studied 
theoretically  and  leads  to  the  conclusion  that  above  Mach  2  the  thickness 
would  deviate  substantially  from  the  Navier-Stokes  predictions.  However,  all 
of  these  theories  suffer  from  convergence  problena  in  the  series  solutions. 
An  alternative  approach  was  proposed  by  Mott.Smith,  namely  to  uses  a  bimodal 
velocity  distribution,  and  it  has  served  as  a  basis  for  a  nuafeer  of  further 
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studies.  The  blnodal  theory  predicts  substantially  thicker  shock  fronts 
than  the  Navler-Stokes  equations,  particularly  in  the  region  above  Mach  2. 

The  situation,  therefore.  Is  that  the  real  differences  between  the 
various  theoretical  approaches  develop  above  Mach  2  where  there  are  very 
few  measurenents.  Consequently,  it  seened  important  to  extend  the  masure> 
isents  of  shock  front  thicknesses  to  shock  strengths  substantially  above 
Mach  2. 

It  appeared  to  us  that  of  the  three  experimental  techniques  currently 
used  for  studying  shock  structure,  the  optical  reflectivity  method*’**  had 
the  best  chance  of  success  at  higher  Mach  nufri>ers.  Although  the  electron 

lO 

beam  technique  has  been  used  to  Investigate  shock  structure  up  to  Mach  6  , 

the  results  obtained  were  only  of  a  preliminary  nature  and  were  difficult 
to  interpret.  As  for  the  free  molecule  probe  method*  *,  the  problems 
involved  in  producing  stationary  shock  fronts  beyond  liach  2  set  an  upper 
limit  to  the  range  of  shock  strengths  which  could  be  studied. 

In  order  to  make  optical  reflectivity  measurements  on  stronger  shocks 
than  heretofore,  it  was  necessary  to  strengthen  the  apparatus  and  to  employ 
such  higher  driver  pressures  than  in  earlier  work,  since  the  requirement 
that  the  thickness  be  substantially  less  than  the  wavelength  of  the  light 
employed  precluded  the  lowering  of  the  initial  pressure.  In  mny  eases 
combustion  driving  was  employed.  We  were  thus  able  to  generate  shock  waves 
as  strong  as  M  -  4.7  with  initial  pressures  of  approximately  one  atmosphere. 
Furthermore,  increased  aperture  and  better  baffling  of  the  optical  system 
significantly  Increased  the  signal»to*noise  ratio,  making  it  possible  to 
extend  the  measurements  to  lower  reflectivities.  In  this  way,  we  were  able 
to  make  measurements  up  to  Mach  4.70  with  greater  precision  than  previously 
achieved. 
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Another  problem  In  which  we  were  interested  was  whether  rotational 
relaxation  in  a  diatomic  gas  was  as  rapid  in  strong  shocks,  where  the 
final  temperature  is  high,  as  it  has  been  found  to  be  in  weak  shocks. 

The  experimental  thicknesses  for  weak  shocks  can  be  accounted  for  satis* 
factorlly  if,  in  addition  to  the  shear  viscosity,  the  gas  also  has  a  bulk 
viscosity.  In  this  bulk  viscosity  is  equivalent  to  a  rotational  relaxa* 
tion  time  of  3*5  collisions.  However,  in  both  Nj  and  CO  the  apparent  final 
density  was  about  51^  lower  than  expected  at  M  =  1.55,  implying  that  in  the 
shock  front  the  "effective”  heat  capacity  was  about  931^  of  the  equilibrium 
value  for  translation  and  rotation.  More  recently,  preliminary  measure* 
ments  of  the  reflectivity  of  shock  fronts  in  oxygen  at  Mach  3^  indicated 
an  even  greater  discrepancy  between  the  observed  and  the  theoretical 
equilibrium  densities.  Since  these  studies  implied  that  not  all  of  the 
rotational  heat  capacity  was  effective  in  the  shock  front  region,  it  was 
of  interest  to  undertake  a  more  extensive  study  of  translational*rota* 
tional  energy  exchange  in  a  diatomic  gas  at  high  temperatures. 

The  Optical  Reflectivity  Method 

In  the  optical  reflectivity  method  a  eollimated  beam  of  light  from 
an  intense  high*pres<,ure  mercury  lamp  is  reflected  off  the  shock  front, 
passes  through  a  narrow  band  interference  filter,  and  is  brought  to  a 
focus  on  the  cathode  of  a  photomultiplier  tube.  The  reflectivity,  R, 
is  given  by  the  equation, 

R  =  R'Wcos  9) 


(1) 


where  An  is  the  change  of  Index  of  refraction  across  the  shock  front* 
9  is  the  angle  of  incidence  (and  reflection)  between  the  light  beam 
and  the  normal  to  the  shock  front*  and  X  is  the  wavelength  of  light. 
R*(X/cos  9)  is  a  reduced  reflectivity  given  by  the  expression* 

R*(X/co8  9)  =  (^n)*  ^  ^  exp[-i<^  cos  9/X)x]dx  |*  (2) 

the  integral  being  the  Fourier  transform  of  (^)*  shock  front 

density  profile  used  in  calculating  R*  in  this  work  was  of  the  form* 


p(x)  =  Pt  +  ^  (3) 

1  +  exp(-4x/L) 

where  L  is  the  maximum  slope  thickness.  This  profile  is  exact  for  the 
bifflodal  theory  of  MottoSmith  and  approximates  all  other  shock  front 
models  closely  for  low  Mach  numbers. 


An  =  (no  -l)(Ji)(^) 
°  Po  Pi 


(4) 


where  the  subscript  "o"  refers  to  N.T.P. *  and  the  index  of  refraction 
varies  with  x*  the  distance  along  the  shock  tube*  in  exactly  the  same 
manner  as  does  the  density.  The  reduced  reflectivity  thus  becomes 


fn^L  cos  0/X 

*  ^sTnKTnaL^os'-STX) 


T 


(5) 


By  varying  the  initial  pressure*  pi*  and  by  using  the  relation 

L  =  Lo(f2)*  where  is  the  thickness  at  one  atmosphere*  one  may  get 
Pi 

an  experimental  curve  for  R*  as  a  function  of  (Xpi/cos  0).  At  low 
values  of  R**  this  curve  is  very  steep  and  enables  one  to  determine  Lq 
with  considerable  accuracy.  At  high  pressures*  R*  approaches  unity  and 
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R  approaches  (An)  (1  -i-  tan  d)/4ni  •  Thus,  the  density  change  Imnedlately 
across  the  shock  front  may  be  directly  measured.  This,  in  turn,  affords  a 
means  of  studying  whether  rotational  relaxation  processes  are  completed 
in  the  shock  front. 


Description  of  Apparatus 

The  shock  tube  was  made  of  2.76  in.  i.d.  by  0.120  in.  thick  stainless* 
steel  tubing  and  consisted  of  a  5  l/2  foot  long  high  pressure  section  and 
11  l/4  foot  long  low  pressure  section.  The  high  pressure  section  was 
sufficiently  long  to  prevent  the  overtaking  of  the  shock  wave  by  the 
reflected  rarefaction  wave  when  efficient  driver  gases  were  used.  Because 
strong  shock  waves  were  needed,  provisions  were  made  in  the  high  pressure 
section  for  combustion  driving.  A  7^  helium,  20^  hydrogen  and  10^ 
oxygen  mixture  was  used  in  all  combustlon*driven  shocks.  The  mixture 
was  Ignited  by  the  simultaneous  firing  of  nine  spark  plugs  which  were 
equally  spaced  along  the  length  of  the  high  pressure  section.  These  were 
pulsed  by  the  discharge  of  a  condenser  bank  through  automotive  spark  colls. 
The  apparatus  was  carefully  designed  to  prevent  the  occurrence  of  deton* 
atlons  and  considerable  attention  was  paid  to  the  strength  requirements 
of  the  chamber.  Initial  pressures  up  to  300p.s.l.a.  in  coiid>ostion* 
driven  shocks,  and  up  to  1200  p.s.l.a.  in  hydrogen*drivcn  shocks,  were 
repeatedly  used  without  any  mishaps.  The  shocks  were  initiated  by  the 
spontaneous  rupture  of  scribed  aluminum  dlaphragam,  separating  the  two 
sections  of  the  tube.  Because  of  the  extreme  care  taken  in  scribing 
the  diaphragms,  reaurkable  reproducibility  in  bursting  pressure  was 
obtained.  The  mean  deviation  in  shock  speed  in  about  80^  of  the  expert* 
ments  was  less  than  0.5^  and  was  almost  never  greater  than  1^. 
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Shock  wave  velocities  were  measured  by  means  of  light  screens.  The 
technique  employed  differed  from  the  more  conventional  one^’  in  that 

12.  J.  Emrich,  Tech.  Kept.  II  -  2NR  06l>021,  Princeton  Univ. (19^){ 

W.  Bleakney,  High  Speed  Aerodynamics  and  Jet  Propulsion,  Vol  IV, 

159*  Princeton  Univ.  ^ress,  Princeton  ^195^)» 

light  was  reflected  off  the  front  of  the  shock  wave  rather  than  the 
back.  The  pulses  from  the  light  screens  were  fed  into  a  Beckman  time* 
interval  meter.  The  combined  timing  error  resulting  from  the  rise  time 
of  the  pulses  and  from  the  accuracy  of  the  timer  was  negligible,  being 
of  the  order  of  a  few  tenths  of  a  percent. 

A  schematic  diagram  of  the  optical  system  employed  in  the  reflect* 
Ivlty  measurements  is  shown  in  Figure  1.  The  light  source  was  an  Osram 
HBO  lOOw/2  high-pressure  mercury  arc  lamp.  The  entrance  beam  was  colli¬ 
mated  by  a  simple  convex  lens  and  its  area  was  defined  by  a  rough  pair 
of  slits.  The  beam  then  passed  through  a  shock  tube  window  which  was 
coated  to  minimize  reflection  losses.  The  reflected  beam  traversed  a 
similar  window  and  then  a  set  of  accurate  slits  which  served  to  define 
the  area  of  the  front  from  which  reflected  light  was  accepted.  For 
perfectly  parallel  light,  this  area  was  9  mn  high  and  5^  mm  wide 
(at  9  =  82"),  the  vertical  edges  on  each  side  being  only  9«5  m"  way 
from  the  corresponding  shock  tube  windows.  The  effective  width  of  the 
beam  was  much  smaller,  however,  since  only  the  central  portion  was  very 
Intense.  The  exit  slits  were  followed  by  a  Baird-Atomic  narrow-band 
Interference  filter  centered  at  either  435  <><4  or  579  oM*  Each  of  the 
filters  had  a  7  bandwidth  and  a  peak  transmission  of  about  701^* 
Transmission  at  other  wavelengths  was  less  than  0.1^  of  the  peak  values. 


After  passing  through  either  of  these  filters,  the  bean  was  focused  on 
the  "AO”  slit.  This  slit  was  an  approximately  circular  aperture  and 
determined  the  angular  dispersion  of  the  light  seen  by  the  931A  photo* 
nultlpller  tube  mounted  Immediately  behind  it. 

All  the  optical  components  were  centered  about  a  horizontal  axis 
parallel  to  the  plane  of  the  shock.  In  order  to  eliminate  errors  arising 
from  linear  attenuation  of  shock  strength  and  to  minimize  errors  due  to 
non-linear  attenuation,  the  optical  axis  was  chosen  so  as  to  make  the 
point  of  maximum  reflection  coincide  with  the  midpoint  of  the  velocity 
interval  for  an  angle  of  82*.  Since  the  light  screen  shock  detection 
method  was  also  based  on  reflection  off  the  front  of  the  wave,  this 
arrangement  consisted  in  effect  of  three  velocity  stations  spaced  at 
equal  intervals*  The  peak  of  the  reflectivity  pulse  always  coincided, 
to  within  the  combined  accuracy  of  the  Beckman  timer  and  oscilloscope 
time  scale,  with  the  midpoint  of  the  larger  velocity  interval.  This 
Indicated  that  no  significant  attenuation  of  shock  strength  occurred 
between  the  two  light  screen  stations. 

The  output  pulse  of  the  photomultiplier  tube  was  taken  through  a 
White  cathode  follower  and  fed  into  a  Tektronix  type  5^3  oscilloscope. 

The  signal- to- noise  ratio  was  usually  more  than  twenty.  The  noise 
arose  almost  entirely  from  extraneous  light  reaching  the  photomultiplier 
tube.  For  Mach  numbers  above  3»  an  additional  "noise"  appeared  as  a 
result  of  the  gas  behind  the  shock  front  becoming  luminous  at  the  high 
temperatures  generated  by  such  strong  shocks  (e.g.  about  i|000*K  at 
M  =  5)*  At  low  initial  pressures,  the  luminous  background  was  sufficient 
to  mask  the  reflectivity  pulse  and  it  became  necessary  to  increase  the 
pressure  in  order  to  obtain  adequate  reflection-to-luminoslty  ratios. 
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The  low  pressure  section  of  the  shock  tube  was  evacuated  to 
100  microns  or  less  before  gases  were  admitted.  The  leak  rate  was 
negligible  in  view  of  the  high  initial  pressures  employed.  The  gases 
used  were  taken  from  coaraercial  cylinders  without  any  further  purifi* 
cation.  Pressures  below  an  atmosphere  were  measured  on  a  mercury  maiio« 
meter  and  above  an  atmosphere  by  a  l6  in.  Heise  0-75  P*s.i.  gage  which 
was  accurate  to  within  0.1^. 

The  Experimental  Procedure  and  Method  of  Calculation 

o 

In  all  the  experiments  reported  here.  X  was  either  433^  A  or 
o 

5790  A.  and  6.  with  one  exception,  was  82*.  In  any  particular  group 
of  experiments,  the  wavelength  and  the  Initial  pressure  were  fixed  at 
predetermined  values  for  all  the  runs,  whereas  average  values  were  used 
for  the  temperature.  T.  the  Mach  nun^er.  M.  and  the  maximum  oscilloscope 
deflection  of  the  photomultiplier  signal.  D.  D  is  related  to  the  exper- 
imental  parameters  by  the  expression 

D  =  kR  =  vO  UV*  ^(An)*R«(  .  (6) 

Hnt»  cos  0 

where  k  is  a  constant  which  depends  on  the  intensity  of  the  incident  beam. 

the  transmission  coefficient  of  the  optical  system,  and  the  sensitivities 

of  the  oscilloscope  and  photomultiplier  circuits. 

The  constant  k  could  be  determined  by  making  use  of  the  fact  that 

when  Xpi/cos  6  is  great  enough.  R*  =  1.  Consequently,  it  was  measured 

by  determining  the  deflections  observed  from  Mach  1.^  shocks  in  argon 

at  an  initial  pressure  of  6o  p.s.i.a.  Under  these  circumstances,  using 

values  of  determined  by  previous  Investigators.  R*  was  0.93  end  0.96 

o  o 

for  wavelengths  of  4350  A  and  5790  A  respectively.  Moreover.  R*  was  so 
high  that  it  was  insensitive  to  errors  in  the  value  of  Lq  assumed. 
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In  early  measurements^  this  value  of  k  was  used  to  calculate  the 
average  absolute  value  of  R*  for  a  set  of  experiments  from  Equation  (6), 
and  Lq  was  then  determined  from  Equation  (5)*  If  R*  vas  much  less  than 
unity  this  method  should  have  yielded  accurate  values  of  Lg,  since  R* 
is  then  very  sensitive  to  Lg.  However,  it  was  found  that  the  apparent 
value  of  the  thickness,  Lg,  differed  for  shocks  at  the  same  Mach  nund^er 
but  different  initial  pressures*  It  was  thus  apparent  that  k  was  not 
a  constant  but  depended  on  the  initial  pressure,  pa.,  the  shock  strength, 
M,  or  both*  To  cope  with  this  problem,  reflectivity  measurements  were 
made  at  a  number  of  values  of(X  p^/cos  6),  and  hence  R',  each  set  being 
calibrated  against  shocks  in  argon  at  M  =  1*4  and  p^  =  6o  p*s*i*a*  The 
resulting  data,  as  can  be  seen  from  the  curves  in  Pig*  2,  could  be  fitted 
to  a  theoretical  curve  of  the  form  of  Equation  (5)  to  within  the  experi* 
mental  error*  However,  instead  of  approaching  the  asyn^tote,  R*  =  1, 
they  approached  a  value  g  <  1*  Hence,  they  acted  as  if  Eq.  (5)  were 
correct  but  the  true  calibration  constant  were  gk,  where  g  is  a  function 
of  the  shock  strength  but  not  of  pt*  The  dependence  of  g  on  the  Mach 
number  is  shown  in  Fig.  3*  problems  this  raises  will  be  discussed 

at  greater  length  later* 

For  several  shock  strengths,  measurements  were  carried  out  In  argon 
at  only  one  value  of  the  parameter  (Xpi/cos  9).  For  these  cases,  g  was 
obtained  by  Interpolation  from  those  experiments  In  which  it  was  deter* 
mined  by  the  aforementioned  procedure  (i*e*  from  Fig*  3)*  calibra¬ 

tion  constant  was  then  taken  as  gk  to  obtain  R*  from  Bq*  v6)  and  Lg 
from  Eq*  (5), 

Results 

A  summary  of  the  pertinent  data  and  results  for  the  individual  sets 
of  experiments  in  argon  and  nitrogen  is  given  in  Table  1*  The  error 
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limits  shown  for  gR*  are  the  90%  confidence  levels  for  the  mean  of  this 
quantity.  With  the  two  exceptions  indicated  in  the  table,  all  the 
reflectivity  measurements  were  made  at  an  angle  of  incidence  of  82*  and 
at  a  horizontal  slit  width  of  6.96  mm. 

Experimental  points  obtained  at  approximately  the  same  shock 
strength,  but  at  different  values  of  the  parameter  are  fitted 

to  the  theoretical  reflectivity  curve  as  shown  in  figure  2.  Tables 
2a  (argon)  and  3  (nitrogen)  give  the  values  of  g  and  determined 
from  these  plots.  The  latter  quantity  is  given  in  the  reduced  form, 
(^),  where  li  =  ^  ^  is  the  Maxwellian  mean-free-path  ahead 

of  the  shock.  It  should  also  be  noted  that  the  g  factors  are  normalized 
relative  to  the  reference  value  at  Mach  1.4  which  is  set  equal  to  unity. 

Whereas  Figure  3  shows  g  is  a  relatively  smooth  function  of  M, 
there  is  no  correlation  with  p^,  of  which  it  might  also  have  been  a 
function,  g  factors  corresponding  to  the  Mach  numbers  at  which  only 

Xpi 

one  value  of  (^os  '6^  investigated  were  obtained  by  interpolation 
from  a  smooth  curve  drawn  through  the  points  on  the  argon  plot. 

Results  for  these  experiments  are  given  in  Table  2b.  Since  the  highest 
shock  strength  at  which  g  was  determined  experimentally  was  Mach  4.72, 
a  linearly  extrapolated  value  of  g  was  used  for  the  point  at  Mach  4.85* 

Reciprocal  shock  thicknesses  for  argon  obtained  in  this  investiga¬ 
tion  are  plotted  against  Mach  number  in  figure  4.  Also  sho%m  are  exper¬ 
imental  points  determined  by  other  workers  as  well  as  the  theoretical 
Navier-Stokes  curve  for  the  density  thickness  and  that  calculated  from 
the  bimodal  theory  using  several  intermoleeular  potentials.  Bach  of 
the  theoretical  corves  uses  potential  parameters  which  duplicate  the 
experimental  viscosities  up  to  the  temperatures  reached  in  a  Mach  3 
shock  wave.  Reciprocal  shock  front  thicknesses  listed  in  Table  2a 


11- 


are  denoted  in  the  figure  by  open  circles  and  are  the  most  reliable 
points  since  they  represent  the  greatest  number  of  experiments.  The 
g  factor  for  each  of  these  points  was  determined  at  the  same  Mach 
number  at  which  the  thickness  was  measured,  whereas  for  the  experiments 
reported  in  Table  2b  the  g  factor  had  to  be  estimated  by  means  of  a 
graphical  Interpolation.  It  can  be  seen  from  figure  4.  however,  that 
the  elimination  of  either  set  of  points  would  not  change  the  essential 
form  of  the  experimental  thickness  curve. 

Finally,  results  for  the  experimental  shock  thicknesses  in  nitrogen 
are  compared  in  figure  5  to  the  Navier-Stokes  density  thicknesses,  with 
and  without  a  bulk  viscosity  term.  Each  of  these  points  is  based  on 
measurements  at  three  values  of  (\pi/cos  9). 

The  g  Factor:  Argon 

The  g  factor  may  be  interpreted  as  the  ratio  of  the  actual  cali¬ 
bration  constant  for  the  experiment  to  that  at  Mach  1.4  in  argon.  From 
figure  3»  it  is  seen  that  this  ratio  differs  from  unity  and  is  a  function 
of  the  Mach  number. 

Such  a  situation  can  occur  when  the  sensitivity  of  the  apparatus 
changes  with  shock  strength.  For  example,  if  the  electronic  rise  time 
were  comparable  with  the  rise  time  of  the  experiment,  the  distortion  of 
the  photomultiplier  signal  would  increase  with  shock  speed.  This  parti¬ 
cular  effect  could  be  neglected  in  our  experiments,  however,  since  the 
combined  rise  time  of  the  circuits  and  of  the  oscilloscope  input  leads 
was  of  the  order  of  O.05  iMec  whereas  the  rise  time  of  the  reflectivity 
pulses  was  always  greater  than  one  Usee. 

The  g  factor  would  also  differ  from  unity  if  a  Mach  number  dependent 
property  of  the  shock  wave  such  as  the  density  ratio  or  the  shape  of  the 
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density  profile  differed  from  that  expected*  In  view  of  the  great 
precision  to  which  the  validity  of  the  Rankine-Hugonlot  equations  has 
been  directly  verified  by  very  many  distinctly  different  experimental 
techniques**’*^®**'’  it  seems  most  unlikely  that  (Ap/pi)  could  differ 
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Curvature  has  Indeed  been  observed  in  shock  tubes*  Lin  and  Fy£e»^* 

15.  S.  C.  Lin  and  W.  I.  Fyfe,  Phys.  Fluids,  4,  238  (I96I). 

using  an  ultraviolet  absorption  technique,  measured  shock  curvatures  in 
air  in  a  24"  shock  tube  at  initial  pressures  between  0.02  mm  Hg  and 
0*3  mm  Hg  with  shocks  ranging  from  Mach  12  to  Mach  22*  Duff  and  Young^* 


16.  R.  E.  Duff  and  J.  L.  Young,  III,  Phys*  Fluids.  4,  8l2  (I961). 


measured  the  curvature  of  Mach  1.8  to  Mach  6.3  shocks  in  argon  at 

initial  pressures  of  0.1  mm  to  20  nm  Hg,  using  a  1  1/8"  shock  tube. 

Both  sets  of  investigators  found  the  curvature  to  be  proportional  to 
«l/2 

Pt**  and  independent  of  the  Mach  number.  Lin  and  Fyfe  also  argued 
that  the  radius  of  curvature  is  proportional  to  r*^/^,  if  r  is  the 
radius  of  the  shock  tube. 

Extrapolation  of  these  results  to  our  conditions  predict  a  radius 
of  curvature  of  about  Z^  meters,  rather  than  the  value  of  about  one  meter 
required.  Furthermore,  neither  the  pi"  1/2 

pressure  dependence  or  the 

independence  on  Mach  number  is  reflected  in  the  behavior  of  the  g  factor. 
All  of  these  factors  cast  doubt  on  our  interpretation  of  its  origin. 

It  is  essential  to  note,  however,  that  the  laeasurements  were  all 
made  at  low  pressures  under  conditions  where  the  boundary  layer  was 
laminar.  At  the  Initial  pressures  employed  in  the  present  work  the 
boundary  layer  is  partially  turbulent,  de  Boer  has  studied  this 


17.  P.  C.  T.  de  Boer,  Tech.  Note  BN.297,  UnW.  of  Md.  C1962). 

situation  theoretically!  he  predicts  that  the  curvature  should  increase 
with  Mach  number  and  be  almost  Independent  of  initial  pressure,  in 
agreement  with  our  observations.  Moreover,  the  magnitude  of  the  axial 


bulge  he  predicts  Is  very  much  greater  than  that  extrapolated  from  the 
low  pressure  data.  Another  point  of  interest  is  his  prediction  that 
non-uniformities  of  the  tube  walls  increase  the  curvature.  He  found 
experimentally  that  extremely  small  irregularities  in  the  shock  tube 
surface  and  in  the  alignment  of  the  observation  windows  with  the  walls 
of  his  rectangular  shock  tube  led  to  an  appreciable  increase  in  shock 
curvature. 

An  even  larger  number  of  irregularities  which  could  contribute  to 
shock  wave  disturbance  were  present  in  the  shock  tube  used  in  our  exper¬ 
iments.  Two  oval-shaped  windows  were  inserted  in  our  cylindrical  shock 
tube  so  as  to  make  the  long  dimension  approximately  flush  with  the  shock 
tube  wall.  As  a  result^  the  internal  diameter  of  the  tube  was  decreased 
by  1.7^  at  the  center  of  the  windows.  Furthermore*  the  small  dimension 
of  the  window  protruded*  at  the  center*  almost  l/l6  in.  past  the  shock 
tube  wall*  a  protrusion  that  did*  however*  becoam  smaller  in  successive 
experiments  as  a  result  of  continuous  erosion  by  passing  shock  waves. 
Still  another  possible  cause  of  non-uniformity  in  the  fluid  flow  was 
the  Indirect  result  of  the  slight  collapse  of  the  shock  tube  wall  during 
hell-arc  welding  of  the  window  bosses.  The  section  was  bored  out  but 
an  unexpected  "chattering"  of  the  boring  tool  produced  some  roughness 
in  the  wall  which  was  diminished*  but  did  not  quite  disappear*  with 
subsequent  samothing.  Thus*  from  both  theoretical  and  experimental 
considerations*  the  suggestion  that  the  g  factor  is  caused  by  shock 
front  curvature  seems  reasonable. 

The  g  factor:  Nitrogen 

Using  the  calibration  constant*  k*  determined  with  Mach  1.4  shocks 
in  argon  and  the  theoretical  values  of  Ap/p^*  for  rotationally  equili¬ 
brated  nitrogen*  the  asymptotic  value  of  R*  also  approaches  a  value 
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which  Is  less  than  one.  In  this  case  the  deviation  might  be  ascribed  to 
a  difference  between  the  experimental  Ap  and  the  equilibrium  value  as  well 
as  to  its  source  (curvature?)  in  argon.  There  is  nothing  in  our  experi* 
ment  jfhich  can  distinguish  the  two.  However*  if  it  is  assumed  that  the 
curvature  is  the  same  in  nitrogen  as  in  argon*  so  that  the  g  factor  from 
this  source  is  the  same  in  both  gases  we  are  led  to  the  following  experi¬ 
mental  values  of  (pa)  exptl: 


Mach  No. 

2.01 

2.54 

2.83 

3.34 

3.72 

(Ps)  exptl 

(ps)  theor. 

0.93 

0.90 

1.01 

0.93 

0.98 

On  the  average  the  densities  are  soam  51^  below  the  value  for  rotational 
equilibrium.  However*  we  can  certainly  conclude  that  (pt)exptl*  «  (pt) 
theor.  to  within  the  coofcined  uncertainty  of  the  curvature  assumption  and 
the  experimental  scatter.  In  that  case  we  can  also  conclude  that  although 
in  the  Mach  3.72  shock  wave  the  nitrogen  is  heated  to  1050*K.*  the  gaa 
comes  to  rotational  equilibrium  within  the  shock  front  transition.  In 
this  case  rotational  relaxation  proceeds  essentially  in  parallel  with 
translational  relaxation*  even  when  the  quantum  number  of  the  most  energetic 
quarter  of  the  molecules  changes  from  J  =  12  to  J  =  22. 

The  Thickness  of  Shock  Fronts  in  Argon  and  Nitrogen 
The  reciprocal  shock  thicknesses  in  argon  and  nitrogen  are  shown  in 
Figures  4  and  5*  ^he  reciprocal  thickness  in  argon  is  seen  to  rise  rapidly 
to  a  maximum  at  about  Mach  3*5  end  to  decrease  slowly  thereafter.  The  aaam 
qualitative  behavior  is  observed  for  nitrogen  although  in  that  caae  the 
shock  fronts  are  thinner  and  the  maximum  is  attained  at  a  higher  shock 
strength. 
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For  comparison  we  include  the  Nawier.Stokes  density  thickness  for 
argon  calculated  by  Schwartz.^**  In  this  calculation  he  has  utilized  a 
viscosity.temperature  relationship  which  fits  experimental  data  for  argon 
up  to  1100*K  (MiO)  and  which  is  theoretically  extrapolated  to  higher  temp, 
eratures  by  utilizing  a  Lennard. Jones  intermolecular  potential.  It  is 
clear  that  while  the  Navier.Stokes  curve  has  the  correct  qualitative 
behavior  it  predicts  substantially  thinner  shocks  than  are  found  exper. 
imentally  when  M  >  2.  Also  shown  are  the  bimodal  theory  curves  of  Muckcnfuss*^* 


18.  C.  Muckenfuss,  Phys.  Fluids,  1325  (19^2). 


for  several  intermolecular  potentials  all  of  which  reproduce  the  experi. 
mental  viscosity.  The  bimodal  curves  are  in  better  agreement  with  the 
experimental  thicknesses  than  the  Navier-Stokes  curve.  The  reflectivity 
results  are  seen  to  agree  best  with  the  exponential. repulsive  (exp.rep) 
model.  The  scatter  in  the  experimental  points  and  the  approximate  nature 
of  both  the  fitting  procedures  and  theoretical  calculations,  make  it  difficult 
to  conclude,  however,  that  our  results  unambiguously  confirm  the  predictions 
of  any  one  particular  Intermolecular  potential. 

In  Figure  5»  the  experimental  nitrogen  thicknesses  are  coaq>ared  with 
theoretical  density  thicknesses  calculated  by  Schwartz^*  from  the  Navler.Stokes 

19.  L.  M,  Schwartz,  private  communication  (1963). 


equations  using  experimental  shear  viscosities  and  bulk  viscosity  set  at  zero 
and  at  2/3  of  the  shear  viscosity.  This  value  of  the  bulk  viscosity  corresponds 
to  a  rotational  relaxation  time  of  5*5  collisions. 
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Since  the  agreement  with  the  Navier>Stokes  equations  when  the  bulk 
viscosity  is  included  is  about  as  good  as  was  obtained  with  argon, 
this  value  of  the  bulk  viscosity  must  be  roughly  correct*  One  cannot 
go  further  at  this  time  than  to  say  that  up  to  Mach  3*7  rotational 
relaxation  seems  to  be  completed  during  the  shock  compression  and 
that  the  effect  on  the  structure  of  the  front  can  be  described 
approximately  by  a  bulk  viscosity  term. 
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TABLE  la 

Data  and  Results  of  Individual  Sets  of  Experiments  in  Argon 


o 


Mach 

No. 

No.  of 
Expts. 

No.  of 
Cal lb. 

X,  A 

Pi 

P«  8*  i  •  8* 

Pi 

^*PqCos  0 

X  10* 

_ 

1.70 

l4 

10 

5790 

59.7 

171.24 

.806  ±  .073 

1.70 

15 

11 

5790 

59.7 

171.24 

.877  ±  .092 

1.71 

15 

12 

28.05 

60.45 

.795  +  .092 

1.70 

10 

6 

11.99 

25.84 

.468  A  .017 

1.92 

10 

7 

9.20 

19.83 

.374  ±  .017 

2.02 

12 

8 

5790 

89.68 

257.24 

.775  +  .040 

2.05 

15 

13 

19.55 

42.13 

.640  i  .058 

2.04 

14 

10 

9.98 

19.35 

.386  ±  .053 

2.02 

11 

7 

7.20 

15.52 

.263  ±  .013 

2.20 

11 

6 

7.20 

15.52 

.292  ±  .017 

2.61 

12 

6 

5790 

33.69 

96.64 

.634  ±  .058 

2.61 

12 

12 

15.30 

32.97 

.578  *  .068 

2.61 

11 

8 

7.20 

15.52 

.317  *  .018 

2.81 

10 

7 

7.20 

15.52 

.338  ±  .036 

2.87 

8 

6 

6.21 

13.38 

.256  *  .022 

3.13 

15 

10 

5790 

24.43 

70.O8 

.638  t  .039 

3.07 

15 

8 

5790 

14.65 

42.02 

.619  ±  .057 

3.14 

14 

10 

14.97 

32.26 

.543  *  .061 

*3.09 

13 

14 

l4.6l 

31.89 

.610  *  .073 

•3.18 

9 

15 

16.11 

19.97 

.389  +  .026 

3.09 

11 

6 

7.20 

15.52 

.335  ±  .020 

3.11 

14 

10 

5.79 

12.48 

.179  *  .014 

3.33 

10 

7 

6.21 

13.38 

.284  t  .019 

3.54 

10 

8 

5790 

23.69 

67.95 

.585  i  .093 

3.60 

10 

8 

13.0 

28.01 

.486  *  .037 

3.56 

12 

11 

12.74 

27.45 

.516  i  .052 

3.52 

8 

17 

7.54 

16.25 

.356  t  .043 

3.53 

10 

6 

7.54 

16.25 

.364  i-  .022 

3.64 

8 

9 

7.54 

16.25 

.396  t  .055 
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TABLE  la  (continued) 

Data  and  Reaulta  of  Individual  Sets  of  Experiments  in  Argon 


Mach 

No. 

No.  of 
Bxpta. 

No.  of 
Callb. 

0 

X,  A 

Pi 

P*  8* i • 8« 

X  Pi 

TT*PjjC08  0 

X  10* 

gR' 

3.72 

8 

17 

6.21 

13.38 

.267  *  .033 

3.84 

13 

12 

6.21 

13.38 

.220  *•  .017 

4.03 

11 

6 

5790 

12.0 

34,42 

.470  t  .041 

4.09 

10 

7 

5790 

7.62 

21.86 

.344  ±  .017 

4.04 

11 

8 

9.80 

21.12 

M  +  .033 

4.07 

13 

9 

7.62 

16.42 

.283  *  .027 

4.05 

10 

6 

6.21 

13.38 

.231  +  .015 

4.05 

12 

14 

5.19 

11.31 

.163  ±  .017 

4.26 

13 

9 

7.'*3 

16.01 

.269  +  .023 

9 

9 

7.35 

15.84 

.283  ±  .038 

4.49 

10 

8 

6.21 

13.38 

.182  f  ,016 

4.68 

12 

9 

5790 

8.91 

25.56 

.367  *  .045 

4.73 

16 

16 

8.91 

19.20 

.316  i-  .039 

4.77 

12 

10 

7.33 

15.80 

.195  ±  .027 

4,67 

10 

6 

5.90 

12.71 

.181  +  ,015 

4.73 

7 

7 

5.^13 

11.70 

.170  ±  .018 

4.85 

7 

10 

6.21 

13.38 

0,168  +  ,036 

1*  Beam  slit  width  -  Z»l6  mm 

2.  0  =  76* 
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TABLE  lb 

Data  and  Results  of  Individual  Sets  of  Experiments  in  Nitrogen 


o 


Mach 

No, 

No.  of 
Expts. 

No.  of 
Calib. 

X,  A 

Pi 

P*  Sa 1 • A* 

jL£i _ 

1T*n  COS  0 
X  10* 

gR* 

2.02 

12 

7 

5790 

59.52 

170.7 

0.607  i  .031 

2.01 

12 

10 

16.32 

35.17 

0.529  t  .042 

2.01 

10 

6 

8.29 

17.86 

0.310  *  .014 

2.54 

15 

7 

5790 

44.78 

128.5 

0.518  i-  .057 

2.53 

15 

7 

21.99 

47.38 

0.490  t  .047 

2.55 

12 

7 

7.62 

16. 42 

0.340  i  .033 

2.81 

12 

7 

5790 

29.78 

85.42 

0.592  ±  .P36 

2.86 

10 

7 

15.49 

33.38 

0.550  t  .030 

2.82 

11 

7 

7.15 

15.40 

0.425  ±  .025 

3.32 

14 

9 

5790 

19.44 

55.76 

0.460  t  .051 

3.34 

13 

8 

13.42 

28.92 

0.489  t  .063 

3.36 

10 

9 

5.79 

12.48 

0.302  i  .034 

3.73 

13 

8 

5790 

14.94 

42.85 

0.677  ±  .050 

3.71 

19 

8 

9.36 

20.16 

0.460  i  .040 

3.71 

15 

9 

6.15 

13.26 

0.352  +  .036 
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OyiBLE  2a 


Suonary  of  Esqperlments  In  Argon  Conducted 
at  Several  Values  of  (X  pi/cos  0) 


No.  of  Sets 


Mach 

Ho. 

of  Expts. 
at  dlf. 

Total  No. 
of  Expts. 

Total  Ho. 
of  Calib. 

$ 

R' 

(X,  Pi,  0) 

1.70 

3 

54 

39 

0.871 

0.537 

2.03 

4 

52 

38 

0.776 

0.497 

0.339 

2.61 

3 

35 

26 

0.668 

0.475 

3.11 

8 

95 

78 

0.684 

0.569 

0.489 

0.262 

3.57 

3 

58 

59 

0.603 

0.617 

4.06 

6 

67 

50 

0.537 

0.527 

0.431 

0.303 

4.72 

5 

57 

50 

0.484 

0.402 

0.373 

0.351 

0.196 

0.240 

0.295 

0.282 

0.346 

0.308 

0.306 


table  2b 


Svmngaiy  of  Experiments  In  Argon  Conducted 
at  Single  Values  of  (X  Pi/cos  9) 


Mach  No.  of  No.  of  li 

Ho.  Expts.  Callb.  g^  R*  L 


1.92 

10 

7 

0.808 

0.463 

0.227 

2.20 

11 

6 

0.758 

0.386 

0.259 

2.81 

10 

7 

0.680 

0.498 

0.308 

2.87 

8 

6 

0.673 

0.381 

0.300 

3.33 

10 

7 

0.625 

0.455 

0.332 

3.72 

8 

17 

0.585 

0.456 

0.333 

3.84 

13 

12 

0.573 

0.585 

0.502 

4.26 

13 

9 

0.528 

0.510 

0.302 

4.42 

9 

9 

0.513 

0.552 

0.328 

4.49 

10 

8 

0.505 

0.360 

0.292 

4.85 

7 

10 

0.470 

0.557 

0.291 

^Estimate  from  Figure  J. 
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TABLE  3 


SvuBnary  of  E]q>erlnent8  in  Nitrogen 


Mach 

No. 

No.  of  Sets 
of  Expts. 

Total  No. 
of  E3q>t8. 

Total  No. 
of  Callb. 

li 

L 

®H2 

^Ar 

2.01 

3 

34 

23 

0.231 

0.626 

0.790 

2.54 

3 

42 

21 

0.327 

0.521 

0.732 

2.83 

3 

33 

21 

0.377 

0.700 

1.052 

3 

37 

26 

0.385 

0.507 

0.811 

3.72 

3 

47 

25 

0.393 

0.551 

0.942 

♦Eetlmate  from  Figure  5- 
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Figure  2 


Heasured  Shock  Front  Reflectivities  for  Argon.  Solid  lines  are 
theoretical  curves. 
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Figure  4 


Variation  of  Reciprocal  Shock  Thickness  (Mean  Free  Paths) 

with  Mach  Huaber  (Argon' 

O  Table  2a) 

'  Present  results 
0  Table  2b) 

Q  Hansen  and  Homig'* 

A  Talbot  and  Sherman® 

1.  Havier-Stokes  theory^*,  Lennard- Jones  6-12  potential 

2.  Bimodal  theory^®;  a.  (exp-rep)  potential 

b.  Lezmard-Jones  6-12  potential 

c.  modified  Buckingham  (exp-6)  potential 
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Figure  5 


Variation  of  Reciprocal  Shock  IMckness  (Mean  Free  Paths) 
with  Mach  Ifumber  (Ritrogen) 

O  present  results;  A  Anderson  and  Homig^;  □  Creene  and  Homlg* 
Theoretical  curves  calculated  from  the  Havier-Stokes  equations^® 
for  hulk  viscosity  zero  and  2/5  of  shear  viscosity. 
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